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A  VOLUMETRIC  CONSTITUTIVE  LAW 
FOR  SNOW  SUBJECTED  TO  LARGE 
STRAINS  AND  STRAIN  RATES 


Robert  L.  Brown 


1.  INTRODUCTION 

Snow  is  a  granular  or  porous  material  comprising  an 
ice-air  mixture.  Because  of  its  high  porosity,  it  is  cap¬ 
able  of  undergoing  large  volumetric  deformations  that 
are  largely  irreversible.  Since  snow  is  so  highly  com¬ 
pressible,  its  density  should  naturally  be  expected  to 
be  a  central  property  which  is  used  to  define  other 
properties  such  as  shear  strength  and  thermal  conduct- 
itivity.  Indeed,  a  close  relationship  between  density 
and  many  of  the  mechanical  and  physical  properties  of 
snow  has  long  been  established,  although  density  alone 
is  not  an  adequate  parameter  for  characterizing  the 
mechanical  properties  of  snow. 

Since  density  (or  porosity)  is  an  important  para¬ 
meter,  the  volumetric  properties  of  snow  will  be  ana¬ 
lyzed  in  this  report.  In  particular,  an  attempt  will  be 
made  to  develop  a  volumetric  constitutive  equation 
relating  material  porosity  to  pressure  loading.  Special 
attention  will  be  given  to  monotonic  loads  generating 
strain  rates  ranging  from  a  low  rate  (I O'5  s'  )  to  rates 
above  10  s'1 . 

Snow  is  idealized  as  a  suspension  of  air  voids  or 
pores  in  an  incompressible  elastic-viscoplastic  material, 
and  the  change  in  its  porosity  is  characterized  in  terms 
of  the  collapse  of  these  air  voids  under  compressive 
loading.  This  is  done  by  considering  the  compressive 
loading  of  a  thick-walled  hollow  sphere  of  polycrystal¬ 
line  ice.  This  model,  while  necessarily  a  rather  simplistic 
one,  does  contain  the  essential  features  of  material 
porosity  and  the  properties  of  the  matrix  material.  It 
also  reflects  the  physical  processes  which  are  necessary 
to  produce  material  flow  to  reduce  pore  size,  so  that  the 
basic  form  of  the  resulting  constitutive  equation  should 
be  correct. 

Some  simplifying  assumptions  made  in  the  analysis 
have  to  be  accounted  for.  In  particular,  the  model  used 


here  assumes  no  interaction  between  pores,  whereas 
medium-density  snow  has  a  high  degree  of  intercon¬ 
nection  of  pores.  This  results  in  a  stress  intensif  ication 
in  the  matrix  material  and  reduces  the  material  stiff¬ 
ness  from  what  the  present  analysis  will  indicate.  To 
take  these  effects  into  account,  a  scale  factor  and  a 
work-hardening  term  must  be  added  to  the  derived 
constitutive  law.  This  will  be  discussed  in  more 
detail  later.  Carroll  and  Holt  (1972)  used  a  similar 
approach  to  analyze  successfully  the  response  of 
porous  aluminum  to  transient  loads.  In  their  work, 
they  assumed  the  material  to  be  an  elastic,  perfectly 
plastic  material,*  whereas  in  the  present  study  a  rate- 
dependent  response  of  the  material  must  be  considered. 

This  report  first  considers  the  constitutive  law  for 
ice.  Since  relatively  high  strain  rates  and  predominantly 
compressive  states  of  stress  are  considered  here,  a 
relatively  simple  form  avoiding  tensorial  arguments 
is  used.  The  report  then  considers  a  detailed  develop¬ 
ment  of  the  constitutive  law  for  snow,  followed  by  an 
evaluation  of  the  results  and  a  discussion  of  some  ap¬ 
plications. 

2.  MATERIAL  REPRESENTATION  OF  ICE 

Previous  studies  by  Dillon  and  Andersland  (1967) 
and  Hawkes  and  Mellor  (1972)  have  indicated  that 
polycrystallinc  ice  under  uniaxial  tension  and  com¬ 
pression  possesses  a  rate-dependent  yield  stress.  Ad¬ 
ditional  experimental  work  by  Haynes  (1976)  has  ex¬ 
tended  the  previous  work  to  strain  rates  in  excess  of 
1 .0  s' 1 .  In  tension,  ice  generally  fractures  at  high  strain 
rates,  and  this  critical  stress  is  not  strongly  rate 

*An  elastic,  perfectly  plastic  material  behaves  as  a  linear  elastic 
material  until  a  yield  condition  is  reached.  Thereafter,  under 
continued  loading,  the  material  deforms  in  a  nonhardening 
manner. 


Figure  1.  Variation  of  axial  yield  stress  with  axial 
strain  rate  for  polycrystalline  ice  / from  Dillon  and 
Anders  land  (1967)1. 

dependent.  However,  under  compression,  the  critical 
stress  shows  a  definite  rate  dependency,  even  at  the 
rates  studied  by  Haynes  (1976).  The  present  study 
primarily  concerns  the  compressive  properties  of  ice 
for  intermediate-to-high  deformation  rates. 

Based  on  the  experimental  results  indicated  in  Figure 
1  (Dillon  and  Andersland  1 967,  Haynes  1976),  the  fol¬ 
lowing  constitutive  relation  is  assumed  for  polycrystal¬ 
line  ice: 


S  =  2Ge  S  «  V 

S=V 


(2.1) 

(2.2) 


where  r  and  y  are,  respectively,  the  shear  stress  and 
shear  strain  rate,  C  is  a  material  constant  and  n  is  a 
scalar.  Equations  2. 1-2.3  were  chosen,  since  they 
acquire  a  straight-line  form  on  Figure  1  and  fit  the 
data  better  than  Glen's  law.  Glen's  law  does  work 
well  at  lower  strain  rates,  however,  and  in  many  cases 
is  mathematically  more  tractable  than  the  form  of  eqs 
2.1  -2.3. 

Consider  now  a  compression  test  on  polycrystalline 
ice.  The  only  nonzero  stress  component  is  oxx,  where 
the  x-direction  is  the  axial  direction  of  the  compression 
test.  The  deviatoric  stress  tensor  is 


S  =  o-  —  (tr  a)  1 
-v  -v  3  'v!  'v. 


(2.5) 


where  tr  implies  the  trace  of  the  tensor,  1  is  the  identity 
tensor,  and^  is  the  Cauchy  stress  tensor?^  therefore 
has  the  components 


Sii  = 


2o, 


0  - 


^xx 

3 


0  - 


XX 

3 


(2.6) 


where 


The  rate  of  deformation  tensor  is 


Y  =  SQ+C\n(AD)  AD  »  1 


(2.3) 


{ l  =  \ 


(2.7) 


50,  C,  and  A  are  material  constants;  5,  e,  and  D  are  the 
principal  difference  values  of,  respectively,  the  deviatoric 
stress  tensor,  deviatoric  strain  tensor,  and  deviatoric 
deformation  rate  tensor;  and  G  is  shear  modulus.  For 
instance,  if  5,  and  S3  are  the  maximum  and  minimum 
principal  values  of  the  deviatoric  stress  tensor,  then 
5  =  S,  -S3.  Similar  definitions  hold  for  D  and  e.  The 
matrix  material,  ice,  is  assumed  to  be  incompressible, 
so  that  the  deviatoric  values  of  the  strain  tensor  and 
the  deformation  rate  tensors  equal  the  total  strain  and 
deformation  rate  tensors. 

The  above  equations  are  essentially  those  of  an 
elastic -viscoplastic  material,  so  that  the  material  behaves 
elastically  until  a  rate-dependent  yield  condition  is 
reached,  whereupon  plastic  deformation  ensues.  These 
equations  are  somewhat  different  from  Glen's  flow  law, 


where  j  is  the  gradient  operator  and  v_is  the  velocity 
vector.  Since  the  material  is  incompressible,  the  diver¬ 
gence  of  the  velocity  Field  is  zero, 

^  *  v  =  0  (2.8) 

and  it  can  be  shown  that^>  has  the  components 


y 


dx 


D, 


0 


\_d^\ 

2  3a 


0 


0 


0 


2  dx 


(2.9) 


Ihe  principal  difference  values,  .S  and  D,  are,  respect¬ 
ively, 


D^. 

2  dx 


(2.10) 


of  snow,  particularly  in  the  very  critical  grain  bonds 
where  massive  localized  deformations  take  place.  Thus, 
the  constitutive  equation  finally  formulated  tor  snow 
may  have  to  be  empirically  adjusted  with  a  work-hard 
ening  term. 

This  problem  has  also  been  encountered  in  the  field 
of  powder  metallurgy  where  a  porous  metal  with  an 
elastic,  perfectly  plastic  matrix  material  shows  con¬ 
siderable  work-hardening  characteristics.  Quite  pos¬ 
sibly  the  large  strains  occurring  in  high  stress  regions 
such  as  grain  bonds  are  of  such  an  order  of  magnitude 
that  the  matrix  material  does  work  harden.  Also, 
quite  possibly  the  actual  description  of  the  deformation 
process  is  incomplete.  At  any  rate,  the  work-hardening 
characteristics  of  porous  materials  such  as  powdered 
aluminum  are  not  predicted  bv  experiments  that  have 
been  run  on  the  solid  matrix  material. 


In  Ihe  results  reported  in  I  igure  I,  strains  were  all 
small,  so  that  the  axial  rate  ol  deformation  component 
is  approximately  equal  to  the  axial  strain  rate 
exx.  f  or  uniaxial  compression,  eq  2.3  becomes 

°xx  =  VC1z7(Mz\J.  (2.11) 

By  adjusting  eq  2.1 1  to  fit  the  compression  data 
represented  in  Figure  I 

S0  =  1x1 06  N  m  : 


3.  DEVELOPMENT  OF  THE  VOLUMETRIC  CON¬ 
STITUTIVE  LAW  FOR  SNOW 

Consider  now  the  deformation  of  a  thick-walled 
hollow  sphere  of  an  incompressible  viscoplastic  material 
with  the  constitutive  relation  given  by  eqs  2. 1-2.3.  The 
initial  internal  and  external  radii  of  the  sphere,  a0  and 
b0,  are  chosen  so  that  the  correct  material  porosity 
results.  In  this  paper,  the  density  ratio  is  defined  as 


P 


C=  I.165x106  Nm'1  (2.12) 

A  =  3.3x  10*5s. 

As  can  be  seen  in  f  igure  1,  this  equation  represents  the 
test  results  quite  well.  There  is  an  appreciable  amount 
of  scatter  indicated  by  the  vertical  lines.  This  is  to  be 
expected,  as  there  may  be  a  good  degree  of  variability 
in  freezing  the  ice,  in  forming  the  specimens,  and  in 
testing.  However,  the  average  results  appear  to  be 
well  approximated  by  the  constitutive  law. 

The  constitutive  law,  eqs  2. 1-2.3,  does  not  involve  a 
work-hardening  term.  Previous  experimental  studies 
on  ice  by  Dillon  and  Andersland  (1962),  among  others, 
have  not  indicated  strong  work-hardening*  character¬ 
istics  for  ice.  However,  there  are  no  published  data  on 
ice  involving  large  strains  in  excess  of  several  hundred 
percent.  Therefore,  it  is  quite  possible  that  at  large 
strains  ice  does  exhibit  work-hardening  properties. 

Such  strains  certainly  do  occur  during  finite  compaction 


where  p  and  pm  are,  respectively ,  the  mass  density  of 
the  porous  material  and  the  matrix  mateiial,  ice. 

Under  an  external  time-dependent  pressure  loading, 
the  internal  and  external  radii  j(t)  and  h(t)  change 
with  time.  The  deformation  is  spherically  symmetric, 
so  the  deformed  coordinates  of  a  generic  point  can  be 
expressed  as: 

f  =  '■('O'  t) 

0  =  (3.2) 

0  =  0Q 

where  r0,  0O,  are  the  undetormed  spherical  coor¬ 
dinate  positions,  and  r,  0,  <t>  denote  the  deformed  coor¬ 
dinate  positions. 

Incompressibility  of  the  matrix  material  requires 
the  (acobian  of  the  deformation  to  equal  unity  ;  i.e., 

r‘  dr  _  , 

'isr 1 


3 


•A  work-hardening  material  is  one  that  stiffens  under  plastk 
deformation. 


(3.3) 


Integration  y  iclds 


Fully  elastic  phase 


c3  rl-B(t) 


(3.4) 


In  (his  phase,  the  strains  are  assumed  to  be  infin¬ 
itesimal.  Iheretore  the  three  strain  components  ate 


where  B  is  an  unknow  n  tunction  ot  time.  I  his  equa¬ 
tion  i.an  be  differentiated  twice.  This  results  in 


f ,  ~  r —  ■  f0  (tt  u  ' 
dr 


(3.9) 


r  =  -d  Si  ♦  JL) 

y  rs  3c-’ 


(3.5) 


where  u  is  the  radial  displacement;  and  it  u  is  small 
compared  with  r,  u  r-r 0  may  be  approximated  by  the 
expression 


for  the  radial  acceleration.  The  acceleration  can  be 
expressed  in  terms  of  an  acceleration  potential  C'(r.f) 


u  -SUL. 


(3.10) 


dr 


3'o 


(3.6) 


B\  doing  this,  the  following  can  be  arrived  at  for  all 
j  ^  r  <  b 


for  small  strains,  ru  d if  1  ers  by  a  small  amount  from  r; 
therefore,  in  this  section  the  distinction  will  be  dropped. 
The  strains  become 


**  *■  ^ 


18c 


4  3  r 


(3.7) 


The  above  results  are  strictly  of  a  kincmatical  nature 
and  depend  only  on  the  constraint  of  material  incom¬ 
pressibility.  In  addition,  material  incompressibility  may 
be  used  to  arrive  at  the  following  relations: 


er  =  -  — 
3r3 


(3.11) 


(3.12) 


3r’ 


Since  the  material  is  assumed  to  be  linearly  elastic, 
the  constitutive  equation  acquires  the  form 


(3.8a) 

5r  -  2GYr  4GB/(dr}) 

(3.13a) 

(h3-u3) 

j  _  Oq  (a-1)  hi  _  4° 

(3.8b) 

S0  =  2  Ge0  =  -2Gfl/(3rJ) 

(3.13b) 

(<*0-1 ) 

(<*<)'>> 

=  2Ge0  =  -2Gfi/(3rJ) 

(3.13c) 

Bit) 


‘'o  (ao-°) 

(a0-M 


(3.8c) 


As  the  external  pressure  P  increases,  the  deformation 
proceeds  in  three  distinct  phases: 

(1 )  an  initial  purely  elastic  phase, 

(2)  an  clastic-plastic  phase  with  an  clastic/plastic 
interface  at  c  =  c,  where  a<c  <  b,  and 

(3)  a  fully  plastic  phase. 

During  the  first  two  phases,  the  strains  are  assumed  to 
be  small,  but  during  the  fully  plastic  phase,  large  strains 
can  be  incurred.  Each  phase  is  now  considered  separate¬ 
ly. 


where  5r,  Sg,  and  er,  eg,  are,  respectively ,  the 
deviatoric  stress  and  deviatoric  strain  components.  Since 
the  material  is  incompressible,  the  volumetric  strain  is 
zero,  and  the  hy  drostatic  pressure  p  cannot  be  evaluated 
with  the  constitutive  equation  alone. 

The  radial  equation  of  motion  is 


do,  -s 
_ L  +  — 


{0,-Oq)  =  pmr 


dr  r 

which  may  he  put  in  the  modified  form 

Jp  =  n  W 
dr  m  dr 


(3.14) 


(3.15) 


for  which  the  boundary  conditions  are: 


o,  =  0  r  =  a 

(3. 

or  =  -P(t)  r  =  b. 

(3. 

It  should  be  noted  here  that  pore  pressure  is  assumed 
negligible.  Under  extremely  large  rates  of  loading  or 
near  low  density  ratios  (say  l<a<1.2),  the  wisdom  of 
this  assumption  may  be  questionable.  This  will  be  dis¬ 
cussed  in  more  detail  later.  Integration  of  eq  3.1  5 
results  in: 

-p  (r,t)  =  pmW, f )+/!(/)  (3.18) 

where  h{r)  is  a  time  function.  Since  or  =  5r  -p, 


This  equation  is  the  same  as  Carroll  and  Holt’s  (1972) 
result  in  their  study  of  porous  aluminum. 

Elastic-plastic  phase 

Equation  3.19  indicates  that  the  maximum  stress 
occurs  at  the  inner  radius.  Yielding  then  initiates  there; 
and,  as  the  pressure  P(t)  continues  to  increase,  a  yield 
surface,  r  =  c(t),  propagates  radially  outward.  The  yield 
condition  given  by  eq  2.3  is  reached  at  a  critical  pres¬ 
sure  Z5,  when  the  principal  deviatoric  stress  difference, 
S'-Sq,  reaches  the  critical  value  Y,  which  itself  is  rate 
dependent. 

Assume  now  P  >  P^  and  that  yielding  has  propagated 
out  to  a  radial  distance  c,  a<c<.b.  For  r<c,  an  inelastic 
stress  state  exists,  and  for  r*c,  the  material  is  elastic. 
When  analyzing  the  outer  elastic  zone,  the  following 
boundary  conditions  apply: 


ar=pm^/t+4Gtf/(3r3).  (3.19) 

Applying  the  boundary  conditions  gives  the  solution 

(3.20) 

3  a2  b 2 

ol  can  then  be  expressed  directly  in  terms  of  P{t)  by 
making  appropriate  use  of  eqs  3.7  and  3.8.  These, 
combined  with  eq  3.20,  yield 

P(t)  =  t2  Q(a,  a,  a)*-—0^  (3.21) 

3a(a-l) 

where 

f2  =Pm  «o/l3  K*')2/3I 

Q  =  -or ( or- 1 ) " 1  / 3  -ar" 1  1  (3.22) 

+  la2((a-1)-4/3-«-4/3|. 

6 

Equation  3.22  describes  the  material  volumetric  response 
to  a  hydrostatic  pressure  loading.  The  term  P(t)  re¬ 
presents  the  actual  pressure  in  the  matrix  material. 

Carroll  and  Holt  (1972)  have  shown  that  the  average 
pressure  in  the  porous  material  fi{t)  is  approximately 

P =  P/a.  (3.23) 

As  long  as  the  material  does  not  yield,  eq  3.21  des¬ 
cribes  the  time-dependent  response  of  the  material. 


or  =  -P(t),  r  =  b 

(3.24) 

°r  =  -pc(t)>r~C 

(3.25) 

SrSg  =  2GB=Y,r  =  c. 

C3 

(3.26) 

At  the  interface  of  the  elastic  zone  and  plastic  zone,  the 
stress  must  satisfy  the  yield  condition.  In  the  yielded 
zone,  eq  2.3  gives 

Sr  =  lY=l(S0+C]nAD) 

(3.27) 

which  remains  valid  as  long  as  the  loading  is  monotonic. 

Consider  first  the  yielded  zone,  where  the  constitu¬ 
tive  law  is  defined  in  terms  of  the  deformation  rate 
tensor.  The  principal  rate  components  arc,  if  vq  and 
vq  are  zero, 

D„  =  bvjbr 

(3.28) 

CD 

il 

(3.29) 

.C 

ii 

(3.30) 

The  shearing  components  arc  zero.  Since  the  deforma¬ 
tion  is  viscous,  a  Eulerian  description  can  be  used  here. 
Incompressibility  requires  that  the  divergence  of  the 
velocity  vector  vanish;  i.e., 
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A%  indicated  by  eg  1-4 the  piessure  response  ot  the 
material  consistent  two  parts:  (a)  a  quasi-static  part 
and  (b)  a  dynamic  part  as  represented  by  the  acceler¬ 
ation  potential  term  pm(^b-^/J)/o. 

I  he  acceleration  term  can  easily  be  shown  to  he 


Simplified  equation 

II  V  is  considerably  smaller  than  G,  it  can  be  shown 
that  the  strains  incurred  (and  therefore  o-o  )  will  be 
quite  small  during  the  first  two  phases  of  the  deform 
ation.  If  such  is  the  case,  fhese  two  first  phases  of 
deformation  can  be  neglected  by  representing  the 
material  as  a  rigid  plastic  material. 
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where  a  and  A  are  expressible  in  terms  ot  «  by  eq  3.8. 

I  be  acceleration  term  becomes  significant  only  at 
very  high  load  rates  which  involve  significant  inertial 
effects,  including  shock  waves,  impact,  etc.  Presently 
we  are  concerned  with  quasi-static  load  situations  so 
that  the  acceleration  terms  may  be  neglected. 

As  indicated  in  section  2,  some  modification  ol  the 
constitutive  law  is  needed  to  reflect  work-hardening 
effects  and  to  provide  some  correction  tor  some  ot  the 
simplifying  assumptions  used.  By  adding  such  a  work- 
baldening  term,  /  e\p(-yta/a0),  and  by  making  use  of 
eq  3.24,  eq  3.43  can  be  put  in  the  following  form: 

r  r  2(S0-CK  m  |  /  exp 

3ft  ft- 1  o(a-1) 

(-Oft/ftp)  (3.4  * 

where  0  and  /  are  material  coeff  icients.  I  his  wotk- 
hardemng  term  is  ol  a  nature  consistent  with  the  teinr 
used  by  St  Lawrence  and  Bradley  (1*174),  who  noted 
that  Gilman  (l‘lh'1)  used  a  similar  term  lor  lithium 
fluoride  crystals  to  account  for  variations  in  the  mobile 
dislocation  velocity.  In  porous  snow,  work  hardening 
may  be  due  to  excessive  amounts  ol  inter gianulai 
motion  that  would  not  be  present  in  solid  polycry  slab 
line  ice.  St.  Lawrence  and  Bradley  (1*174)  addressed 
this  question  in  a  more  thorough  manner.  In  addition, 
the  model  used  here  is  rather  ideal  in  the  sense  that 
isolation  ot  the  pm  ex  is  assumed,  so  that  the  effect  ol 
one  pore  on  another  was  neglected. 


♦Clnlb^Li  /  exp(-0iL)  >£  ( 

<»(«- 1 )  o0  1 

where  is  the  critical  stress  at  which  the  matrix 
material  first  becomes  fully  yielded.  is  found  by 
setting!  />  in  eq  3.42  to  yield 


ft  _Lln(-22J|2(S0-O 

1  In.,  «o-’ 
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4.  COMPARISON  WITH  EXPERIMENTAL  DATA 

.V0,  C,  and  4  are  given  in  eq  2  12.  1  he  only  other 
mater  ial  constants  that  need  to  be  evaluated  aie  /  and 
0.  the  two  work  hardening  coefficients.  Based  on  data 
supplied  by  Abele  and  Gow  (1*>76),  these  two  coef¬ 
ficients  were  found  to  be 

/  3.07 


I  q  3.48  was  adtusted  lo  data  (Abele  and  Gow  l*>7b) 
tor  snow  at  -10‘Y  and  with  an  initial  density  of  400 
kg  m  V 

I  he  variation  ol  the  pressure  p  with  ft  can  then  be 
calculated  by  means  of  eqs  2.12,  3.47  and  4.1  and 
compared  with  the  data  supplied  by  Abele  and  Gow 
( 1  *>75,  I97h).  Hie  details  of  the  experimental  work 
are  given  in  the  reports  ot  Abele  and  Gow,  and  the 
reader  is  referred  to  them  tor  a  mote  involved  des¬ 
cription  of  the  experimental  setup.  Briefly  ,  the  tests 
consisted  of  high-rate  axial  compression  tests  in  which 
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the  specimens  were  constrained  laterally  in  Teflon- 
lined  aluminum  cylinders.  As  a  consequence,  the 
deformation  resulted  in  large,  easily  measured  volu¬ 
metric  deformations. 

Consider  now  the  problem  of  the  uniaxial  com¬ 
pression  of  snow.  If  we  take  the  x,  direction  as  the 
axial  direction  and  the  x2  and  x3  directions  as  the  two 
lateral  directions  of  a  cylindrical  specimen,  the  motion 
can  be  described  by  the  equations  (see  Fig.  2): 

x,  =X'J+Xf(1-X',//t0) 

*2  ~x2 

*3=  *3  (4-2) 

where  A", ,  and  are  the  undeformed  coordinate 
positions,  x  ( ,  x2  and  x3  are  the  deformed  coordinate 
positions,  t  is  time,  X  is  the  crosshead  velocity  (constant 
in  this  case),  and  h0  is  the  initial  specimen  height.  The 
Jacobian  of  the  deformation  gives 


y  =£°=  HV*0)f 

p 

and  therefore  the  density  is 

p-poIMWoM'1. 

The  density  ratio,  however,  is 
a  =  a0  ( 1-(X//t0)f| 
and  its  rate  is 

a  =  -^L. 

h0 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


The  constant  rate  tests  do  give  a  constant  density 
ratio  rate.  Note,  however,  that  the  rate  of  change  of 
density  is  not  constant.  For  this  set  of  experiments, 


Figure  2.  Description  of  confined  compressive  tests. 


a  was  readily  calculated  by  eq  4.6,  and  eq  3.47  was 
then  utilized  to  find  the  variation  of  p  and  a  for  a 
given  rate.  Figure  3  shows  the  comparison  of  theory 
with  experiment  for  snow  at  -10°C  and  for  three 
initial  densities.  The  data  acquired  by  Abele  and  Gow 
(1975,  1976)  included  a  variety  of  strain  rates,  since 
different  specimen  sizes  and  crosshcad  speeds  were 
utilized.  The  data  shown  here  reflect  rates  in  the 
neighborhood  of  a  %  10  s'1 . 

The  data  available  from  Abele  and  Gow  (1975,  1976) 
measured  only  o( ,  the  major  principal  stress,  whereas 
the  hydrostatic  pressure  /?=  1/3  (0|+o2+o3)  was 
needed.  Therefore,  a  series  of  experiments  was  run 
to  measure  lateral  stress  as  well  as  axial  stress  so  that 
Abelc’s  and  Gow’s  data  could  be  adjusted  to  reflect 
prather  than  o, .  It  was  found  that 

0.6o,<£n0.98o,  (4.7) 

for  all  the  experiments  run,  thereby  fairly  well  bracket¬ 
ing  the  data  report  by  Abele  and  Gow  (1975,  1976). 

One  interesting  result  was  observed  in  the  testing 
program.  In  tests  involving  "old"  snow,  the  lateral 
stress  was  found  to  be  about  90%  of  the  axial  stress, 
but  for  tests  of  "new"  snow,  the  lateral  stress  was 
only  3040%  of  the  axial  stress.  In  each  case,  unbonded 
snow  (sifted  within  two  hours  of  the  time  of  testing) 
was  used,  so  that  the  difference  between  the  lateral 
stress  in  the  old  snow  and  that  in  the  new  snow  must 
have  been  due  to  differences  in  crystal  structure, 
primarily  crystal  shape.  The  old  snow  had  undergone 
cquitemperature  metamorphism,  and  the  crystals  had 
a  spherical  shape.  The  grains,  therefore,  were  capable 
of  rolling  and  sliding  relative  to  each  other,  thereby 
accommodating  a  lateral  motion  to  produce  a  large 
lateral  stress.  Flowevcr,  new  snow,  with  its  complicated 
grain  structure,  would  not  allow  this,  thereby  producing 
a  smaller  lateral  stress. 

As  can  be  observed  from  Figure  3,  the  comparison 
between  theory  and  experiment  is  quite  good,  with 
the  essential  characteristics  of  the  deformation  being 
represented  by  the  theory. 
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f  igure  3.  Comparison  of  theory  with  experimental  data 
for  snow  at  -IO°C  at  three  initial  densities. 


Figure  4  compares  theoretical  pressure  curves  with 
data  collected  and  summarized  by  Mellor  (1974).  This 
figure  contains  results  of  laboratory  studies  as  well  as 
field  data  relating  density  to  gravity  induced  pressure. 
Once  again  the  theory  looks  quite  reasonable  compared 
with  the  data.  It  should  be  remembered  that  some  of 
the  data  shown  in  this  figure  represent  uniaxial  stress 
conditions,  and  that  the  actual  hydrostatic  pressures 
arc  only  one-third  the  values  shown  for  these  data. 
Therefore,  some  of  the  experimental  curves  would 
move  down  vertically  relative  to  the  theoretical  curves. 
However,  since  there  is  such  a  diverse  range  of  load 
histories,  temperatures,  and  time  ranges  contained  in 
Figure  4,  any  meticulous  adjustments  would  not 
necessarily  change  things  that  much.  What  can  be  said, 
though,  is  that  cq  3.48  appears  to  be  functionally 
correct  for  snow  with  initial  densities  exceeding  300 
kg  m‘3. 

Figures  5  and  6  further  describe  some  important 
properties.  The  deformation  rate  dependency  is  il¬ 
lustrated  in  f  igurc  S.  For  snow  with  an  initial  density 
of  350  kg  m'3,  the  stress  response  is  shown  as  a  function 
of  density-ratio  rate  at  three  different  instantaneous 
densities.  As  can  be  seen,  a  rate  dependency  does  exist, 
but  the  importance  of  rate  decreases  as  rate  increases. 


Figure  4.  Comparison  of  theory  with  laboratory 
and  field  data  for  A  and  B:  Natural  densif n  a¬ 
tion  of  snow  deposits.  F,  |,  and  K.-  Uniaxial 
compression  of  snow.  E.-  Calculated  values 
of  plane  wave  impact. 

At  rates  characteristic  of  stress  waves,  therefore,  one 
might  be  able  to  assume  a  constant  yield  stress  and 
achieve  a  simplified  version  of  the  constitutive  law 
given  in  eq  3.48.  However,  for  lower  rates,  say  in  the 
range  10's  *.  I  ok  1 0,  the  rate  dependency  is  significant 
enough  that  such  a  simplifying  procedure  would  not 
be  recommended. 

Figure  6  gives  an  indication  of  just  how  effectively 
snow  can  absorb  energy  during  compaction.  In  par 
ticular,  this  figure  shows  the  work  required  to  compact 
snow  to  a  terminal  density  of  700  kg  m3  for  a  range 
of  initial  densities  and  density-ratio  rates.  One  can  see 
immediately  that  initial  density  has  a  dramatic  effect 
on  work  required  to  compress  snow.  Density-ratio 
rate  is  also  significant  but  certainly  much  less  so  than 
initial  density. 


5.  APPLICATION  TO  VEHICLE 
MOBILITY  PROBLEMS 

One  obvious  application  of  this  volumetric  con¬ 
stitutive  law  involves  vehicle  mobility  in  snow.  Consider, 
for  example,  a  vehicle  travelling  at  a  given  speed  along 
a  straight  path  in  a  shallow  snowpack.  The  vehicle  is 
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Figure  5.  Variation  of  pressure  response  with  rate  of 
change  of  density  ratio. 


f  igure  6.  Effect  of  initial  density  on  work  required 
to  compress  snow  to  a  density  of  700  kg  m'\ 
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Figure  7.  Profile  and  pressure  bulb  cross-section  descriptions  of  deformation  pro¬ 
duced  by  tracked  vehicle  in  snow.  P*  =  Nominal  track  pressure. 


assumed  to  be  tracked.  By  shallow  snowpack  is  meant 
a  snowpack  that  is  shallow  enough  so  that  the  pressure 
bulb  cross  section  reaches  to  the  ground.  Figure  7 
shows  a  typical  pressure  bulb  cross  section.  Harrison* 
has  observed  that  pressure  bulbs  usually  have  a  cross 
section  that  is  nearly  rectangular  in  shape.  Some  later¬ 
al  spreading  of  the  bulb  does  occur,  but  the  added 
cross-sectional  area  is  in  most  cases  insignificant. 

In  situations  where  the  vehicle  is  travelling  at  a 
uniform  velocity  and  is  not  experiencing  any  significant 
amount  of  slipping,  it  can  be  reasonably  assumed  that 
the  vast  majority  of  the  energy  expended  in  traversing 
a  snow  field  is  associated  with  the  energy  of  compacting 

•Private  Consultation,  W.L.  Harrison,  CRREL. 


the  snow  beneath  the  vehicle  tracks.  Consider  a  vehicle 
that  has  an  idealised  track  loading,  as  depicted  in  I  igure 
7.  The  track  pressure  is  assumed  to  build  up  linearly 
during  the  entry  portion  of  the  track  loading.  Here 
this  is  assumed  to  be  0.5  m,  with  the  rest  of  the  track 
having  a  length  of  5  nr.  The  entry  length  of  0.5  m  is 
somewhat  arbitrary,  since  this  depends  on  the  depth 
that  the  vehicle  sinks  into  the  snow. 

The  vehicle  power  associated  with  snow  compaction 
is  now  calculated.  Since  the  volumetric  strains  generated 
in  the  snowpack  are  large,  such  an  analysis  must  use  a 
constitutive  equation  valid  for  high  strain  rates  and 
large  strains.  Also,  any  additional  power  requirements 
associated  with  vehicle  slippage  and  drag  would  have 
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to  be  added  to  the  lesults  which  follow.  The  intent 
here  is  to  study  the  manner  in  which  snow  absorbs 
energy  through  compaction. 

I  he  stress  power  is  given  by 


f(l.a) 


1 1 (/  |n(-®J)-2(50-O| . 

C  a- 1 


(S.8) 


*  (5.1) 

where  £and arc,  respectively,  the  Cauchy  stress  ten¬ 
sor  and  the  rate  of  deformation  tensor.  //-(-)  is  the  trace 
ot  the  tensor  quantity  inside  the  parenthesis.  The  stress 
power  M’  is  simply  the  rate,  per  unit  mass,  at  which 
work,  is  being  done  internally  by  the  stress  tensor.  I  or 
a  purely  viscous  material  under  isothermal  conditions, 

W’  would  reduce  to  the  rate  ol  energy  dissipation.  I  or 
a  viscoelastic  material,  the  stress  power  would  con¬ 
it  ibute  to  berth  the  t  ate  of  change  of  sti  am  energy  and 
the  rate  ot  energy  dissipation. 

In  a  unidirectional  deformation,  the  stress  power 
associated  with  compaction  is: 


p  3.x 


where  v  is  the  vertical  particle  velocity,  and  x  is  the 
deformed  coordinate  position  of  the  particle.  The 
continuity  equation  tor  unidirectional  motion  is 

^  -  0.  (5.3) 

at  cw 

Equations  5.2  and  5.3  can  then  be  used  to  obtain 

W  =  -- L^a.  (5.4) 

Pm 


Integration  results  in 
t  A 

H '  -f£adt  (5.5) 

'oPm 

In  the  above,  the  pressure  p  is  a  function  of  time, 
f  quation  3.48  may  be  inverted  to  yield 

a  =  a)  (5.6) 

where 

<2(/.a)  a)+ln|a(a-l)|  |-In4  (5.7) 


The  pressure  loading  in  the  material  below  the  track 
is  assumed  to  have  the  form 

P(t)  P‘  (5.9) 

where  //(f)  is  the  Heaviside  step  function.  I  igure  7 
shows  the  nature  of  the  pressure  distribution  under  the 
track.  The  total  work  done  by  the  track  to  a  unit  mass 
of  snow  is  then 

H  ./'•  i-ft, ),*»«.«>,/,  (5.10) 

0  Pm 

where  f  *  is  the  duration  of  time  that  the  snow  is  undei 
the  track,  f.quations  5.7  to  5.10  can  then  be  used  to 
study  the  energetics  of  oversnow  vehicle  travel. 

By  substituting  the  assumed  pressure  function  given 
b\  eq  5.9  into  eq  S.10  and  integrating  over  the  time 
t*.  the  actual  work  in  compressing  a  unit  mass  of  un- 
deformed  snow  is  found.  1  his  gives  a  direct  measure 
of  track  efficiency  for  a  given  set  ot  parameters  such 
as  vehicle  speed,  track  pressure,  track  geometry,  and 
snow  properties. 

Equation  5.2  gives  the  instantaneous  power/unit 
mass  of  snow  while  the  material  is  under  the  track,  and 
this  expression  can  vary  considerably  during  the  inter¬ 
val  f*  of  track  loading.  However,  tor  the  purpose  ot 
this  study ,  the  average  power  is  more  meaningful,  and 
this  is  calculated  by  dividing  the  total  work  pet  unit 
mass  M'  by  the  time  period  t*  requited  to  produce  this 
work.  Therefore, 

r  *  *  »  ,,$U. «l 

rLi  - dt  (5.ID 

f*b  o 

gives  the  power  per  unit  mass  and  is  subsequently  re¬ 
ferred  to  as  "specific  power."  lhe  specific  power  P 
would  have  to  be  augmented  by  a  factor  equal  to  the 
total  mass  of  snow  under  the  tracks,  it  one  wanted  to 
find  the  total  power  requirements. 

Nominally,  a  track  length  of  5  m  was  chosen,  f  rack 
pressure  was  limited  to  5\I04N  m  *.  which  admittedly 
is  uncharacteristically  low  when  considering  military 
vehicles.  However,  such  pressures  are  more  realistic 
tor  vehicles  such  as  snowmobiles.  Initial  snow  densities 
studied  range  front  300  kg  m  '  to  700  kg  m  \  and 
vehicle  speeds  are  restru  ted  to  about  15ms'1, 
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Figure  8.  Variation  of  specific  power  with  vehicle  sfteed  for  track  pressure. 


Figures  8a  and  b  compare  specific  power  require¬ 
ments  for  two  track  loadings.  A  significant  variation 
in  power  results  when  the  track  loading  is  increased 
from  1x1 04  N  m'2  to  5x104  N  m'v  In  particular,  the 
efficiency  of  the  lower  track  pressure  becomes  increas¬ 
ingly  apparent  at  higher  vehicle  speeds.  This  is  due  to 
the  decreasing  amount  of  snow  compaction  that  occurs 
while  the  snow  is  under  the  track  as  the  speed  increases. 
At  the  lower  track  loading,  there  appears  to  exist  a  criti¬ 
cal  spec  i  above  which  little  increase  in  efficiency  is 
achieved  with  higher  speeds.  For  the  higher  pressure, 
signific  tnt  amounts  of  compaction  continue  to  occur 
at  the  higher  vehicle  speeds.  These  high  pressures,  even 
at  high  speeds,  force  the  vehicle  a  significant  distance 
down  into  the  snowpack,  thereby  expending  much 
energy  in  snow  compaction.  At  lower  track  pressures, 
the  compaction  becomes  much  less  significant,  and  the 
vehicle  tends  to  ride  up  on  the  pack,  which  is  a  much 
more  efficient  configuration. 


Notice  the  tremendous  increase  in  power  when  the 
track  pressure  is  increased  from  1x1  (r  N  m‘2  to 
5xl04  N  m'\  Again,  this  can  he  attributed  to  the 
highly  nonlinear  relationship  between  p,  a.  and  a  in  the 
constitutive  equation. 

Figures  ‘)a  and  h  illustrate  the  variation  of  vehicle 
power  with  initial  snow  density.  In  I  igurc  9a,  the 
relationship  between  specific  power  and  density  is  il¬ 
lustrated  for  iwo  different  vehicle  velocities.  In  I  igurc 
9b,  the  same  relationship  is  shown  lor  three  different 
velocities.  As  can  readily  be  seen,  tor  snow  with  initial 
densities  above  TOO  kg  m  ' ,  a  track  loading  of  lxl O’* 

N  m'2  operates  fairly  efficiently.  However,  the  same 
cannot  he  said  for  a  track  loading  of  5x104  N  m  *. 
where  good  efficiencies  are  not  achieved  tor  initial  snow 
densities  below  500  kg  m 


* 

a 

* 


Track  PraMura  5a  lO-  N  »  1 


Iniiioi  Daftti'y  (kg  *r»> 


j.  /rack  pressure  1*1(7*  ,V  m  *. 
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6  CONCLUSIONS 

The  volumetric  constitutive  equation  given  by  eq 
3.48  has  been  shown  to  correctly  represent  the  response 
ot  snow  to  large  monotonic  compressive  deformation. 
This  equation  was  based  on  a  model  of  pore  collapse 
ot  a  hollow  sphere  ot  ice,  thereby  giving  a  dynamic 
pore  collapse  relationship  relating  pressure  to  material 
Porosity.  Since  work  hardening  was  not  included  in 
the  constitutive  equation  of  the  matrix  material  (ice), 
the  equation  had  to  be  altered  to  reflect  the  work- 
hardening  characteristics  of  snow. 

Ihe  results  do  verify  that  the  hydrostatic  pressure 
is  very  dependent  on  initial  porosity  (or  density)  and 
only  slightly  dependent  on  rate,  so  that  a  change  of  tour 
to  five  orders  of  magnitude  in  strain  rate  is  required 
to  produce  an  order  ot  magnitude  change  in  pressure. 
However,  a  small  change  in  initial  density  can  have  a 
very  significant  effect  on  the  pressure  response  of  snow. 

The  valid  range  of  porosity  rates  is  apparently  -10  s 
s  *a<-IO  s'1  tor  the  flow  law  used  here  for  ice.  The 
apparent  valid  range  tor  initial  densities  includes  300  kg 
m  «.p0vbOO  kg  m'J ,  but  at  very  low  densities  (less 
than  300  kg  m'),  it  is  doubtful  that  the  physical  model 
used  to  develop  the  constitutive  law  would  be  valid,  due 
to  a  high  degree  of  pore  interconnectivity. 

The  results  of  this  study  indicate  that  the  volumetric 
constitutive  equation  can  be  used  as  an  engineering 
tool,  since  it  does  appear  to  give  an  accurate  description 
of  the  response  of  snow  to  volumetric  loading  and  is 
not  unduly  complicated  in  its  simplified  form.  But 
some  additional  experimental  work  would  be  needed 
to  adjust  the  equation  to  various  snow  types.  However, 
in  view  of  the  data  already  acquired  by  Abele  and  Gow 
(1975,  1976),  the  amount  ot  additional  data  required 
to  do  this  would  not  be  very  extensive. 
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